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Abstract 

The scattering of 77r — > tttt is studied using the axial anomaly, elastic unitar- 
ity, analyticity and crossing symmetry. Using the technique to derive the Roy's 
equation, an integral equation for the P-wave amplitude is obtained in terms of the 
strong P-wave pion pion phase shifts. Its solution is obtained numerically by an it- 
eration procedure using the starting point as the solution of the integral equation of 
the Muskelshsvilli-Omnes type. It is, however, ambiguous and depends sensitively 
on the second derivative of the P-wave amplitude at s = which cannot directly 
be measured. 



Key words: Dispersion relations; Integral Equations; Chiral anomaly; Unitarity; Chiral 
perturbation theory; Vector meson dominance; Pion form factor; Multiple scatterings 
PACS: 11.30.Rd; 13.40.-f; 11.55.-m; 11.55Fv; 11.55.Hx; 13.60.Le 

*unite propre 014 du CNRS 







1 Introduction 



One of the fundamental calculation in particle theory is the vr*^ — > 77 decay rate [||]. It is a 
combination of Partial Conserved Axial Current (PCAC) and the short distance behavior 
of Quantum Chromodynamics (QCD): 

^(^0 _ ^ ^F,,e^''"-elKe:K (1) 

with 

^77 = TTTF = 0.025Gey-^ (2) 

where e is the electric charge, f^^ = 0.0924 GeV and A^^^ = 3 is the number of color in 
QCD. This calculation is valid in a world where the vr" is massless. Some corrections have 
to be made in order to take into account of the finite value of the pion mass. It turns 
out that the massless pion anomaly formula is in a very good agreement with the pion 
hfe time data p[, implying that the correction due to the physical pion mass in Eq. (0) 
is very small. 

Another Axial Anomaly result is the process 'jn tttt or its analytical continuation 
7 ^ Svr 0. This last process requires more corrections because, for practical consider- 
ation, measurements are done at an energy far from the chiral limit where the anomaly 
formula is applicable. Furthermore the analytical continuation from one process to the 
other is a delicate procedure due to the presence of the complex singularity which is ab- 
sent in the former reaction. The calculation of the process 77r —>■ tth is in itself interesting 
because of future experiments being proposed at various accelerator facilities and also of 
its important role in the calculation of 7r° — >■ 77* [^]. 

The 77r —* tttt amplitude is given as: 

A(7(fc)7r°(po) ^ 7r+(pi)7r-(p2)) = te"""^ e^pouPi.P2rGis,t,u) (3) 

where s, t, u are the kinematical variables for this process and will be defined below. 

In the chiral limit (the zero limit of the pion 4- momenta), the matrix element is given 
by the anomaly equation: 

G(0, 0, 0) = A = = 9.70GeV-^ (4) 

where the zero in the argument of G{0, 0, 0) refers to the chiral limit of the massless pions; 
the number of colors Nc is equal to 3. 

Experimentally, A is measured at an average photon pion energy of 0.4 GeV and 
assuming that there is no momentum dependence in G(s, t, u), it is equal to [^: 

^expt ^ ^2.9 ± 0.9 ± O.bGeV-'^ (5) 

The agreement between experiment and theory is not good but certainly corrections 
will have to be made because measurements made in this experiment are far from the 
chiral symmetry limit. 
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The calculations of this process are usually done within the Vector Meson Dominance 
models (VMD) ^ [^.Recently this process is discussed within the framework of Chiral 
Perturbation Theory (ChPT) to one loop M and also a combination of ChPT and VMD 



10] and the unitarization of the ChPT two loop amplitude |11 



The purpose of this paper is to investigate the scattering of •yn^ tt+tt^ using disper- 
sion relation, elastic unitarity and the knowledge of the P-wave pion pion phase shits. An 
integral equation is obtained and is similar to the Muskelishvilli-Omnes (M-0) integral 



equation used in the pion form factor calculation |12]. The integral equation obtained here 



is, however, much more complicated due to the symmetry of the problem. Its solution 
can only be obtained by a numerical method. 

The pion form factor calculation, using the M-0 integral equation approach, yields a 
pion radius too low by 10% and a modulus of the pion form factor at the p resonance also 
too low by 15%. This is due to the assumption of the elastic unitarity relation which is 
only valid in the low energy region below 1.2 GeV but cannot be true at and above the 
p'{ 1.5 GeV) region. In order to remove this discrepancy one has to use also, as input, the 
pion r.m.s radius and hence, one has either to make an extra subtraction in the dispersion 
relation or to make use of the polynomial ambiguity of the solution of the M-0 equation. 
One would then get not only the correct value of the absolute value of the pion form factor 
at the p mass (i.e. the p leptonic width), but also a complete agreement with the pion 
form factor below 1 GeV The phases of the form factor are of course the experimental 
P-wave TTTT phase shifts due to the solution of the M-0 equation. 

We also face the same problem in the calculation of the scattering 771" — > vr+Tr" . The 
problem could, however, be more serious here than in the pion form factor calculation due 
to the existence of the t and u channels . Not only the first derivative of the dominant 
P-wave amplitude vanishes at the energy squared s = m^, but its second derivative at 
this energy is not accessible to experiments because of the lack of experimental precision. 

A further comphcation is due to the ambiguity of the solution of the integral equation 
obtained here due to the symmetry of the problem. It is related to but is not purely 
the polynomial type. For this reason we cannot make a comparable prediction of the 
77?*^ — > vr+TT" cross section at the p mass or r{p 717) width. The measurement of 
this width could then be used to make a prediction of the energy dependence of the 
771*^ —>■ TT+TT^ away from the p mass and in particular in the low energy region where the 
the first measurement of A, Eq. (^), was made. 



2 Kinematics and Partial Wave Projection 

The kinematics of this process are defined as: s = {k + Po)"^, t = {pi — po)^ and u = 
(P2 — Po)^- Because all particles involved are on shell, one has s + t + u = 3m^. In the 
center of mass system, in terms of the scattering angle 6, we have: 

SttT'^ s 1 / 

t = ^ h -{s-mD^Jl- Ami/ s cose 
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u = -lis-ml)^l- 4ml /s cos 9 (6) 



The partial wave expansion for G{s,t,u) is given as follows |TT|, [13 



G{s,t,u) = Y.Gi{s)Pl{cose) (7) 

oddl 

where 6 is the scattering angle and P/ is the first derivative of the Legendre polynomial. 
Hence the lowest partial wave is: 

3 



G^{s) = — f dnsm^9G{s,t,u) 



In terms of the function G3^(s,t,M) the differential cross section for the process 77?° 



TT + TT is 



[s-ml)^ -JL^sin" 9 \G{s,t,u)\^ (9) 



dcos^ 10247r' sV2 



3 Vector Meson Dominance and Pion Form Factor 

Because our integral equation solution is a more sophisticated and precised approach 
to the vector meson dominance model (VMD) [|14|, where unitarity and dispersion 
relation are extensively used, it cannot avoid the same problems which are presented 
in these models. Namely the solution can only uniquely obtained when the asymptotic 
behavior of the solution is specified. It is then useful to review briefly the VMD models 
and calculations of the pion form factor with and without introducing the contact term. 



3.1 Vector Meson Dominance Models for jtt^ tt+tt process 

Let us consider the VMD models without and with the contact terms for the 77?° tt+tt" 
process as previouly discussed in the litterature p, 0. Without the contact term, the VMD 
model for 77r° — tt+tt" amplitude is: 

\ 2 2 2 

,1 , X A , m,, m. m,, , , , 

G^-%s,t,u) = -{-^ + + -r^) (10) 

3 — s — t m^ — u 

where are the invariant kinematics and A is defined by Eq. (^). With a contact 

term, it can be written as: 

G-'\s, t, u) = + ^ + ^ - 0} (11) 

3 — c m^ — s mp — t m^ — u 

where c is proportional to the strength of the contact term. Eq. ( pUD can be rearranged 
to give: 

G-'-{s,t,u) = + + [s^t] + [s^ u]] (12) 

3 mi — s 3 — c mi 
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In Eq. (p!0|), the 77r° vr+vr" amplitude vanishes as s,t,u oo while those in Eqs. 
(|Tl|, [l^) do not vanish due to the presence of the contact term. We have introduced 
phenomenologically the contact term c in the scattering 771'' — >■ 7r"'"7r~ without worrying 
how it influences the VMD for the corresponding process P —>■ '-ff. Assuming a complete 
VMD for P — > 77, one has c=l in order that the KSRF relation |]T^ remains valid |^. 



The strength of the contact term also influences the value of the second derivative of 
the P-wave amplitude at s = m^. Expand Eq.(]TTD in a power series of s,t,u: 



ml 1 + + 



G^"^%s, t,u) = X[l + - 1 + + ...] (13) 

6 — c mp 6 — c m^ 



and the P-wave projection of this equation is: 



where the pion mass is introduced by hand. One has flnally: 

d^GAs) , 12 1 



ds^ """"^ 5(3 - c) m4 



A (15) 



Instead of characterizing the contact term by the inflnite energy behavior of the matrix 
element, we can specify its presence by evaluating its second derivative for the P-wave 
at s = m^. For the pure VMD, c = 0, it is equal to (12/15)A(m~^) and for the hidden 
symmetry model 0, c = 1, it is (6/5) A (m 



-4\ 



p 



Eq. ( p!oD yields a decay width T{p iT'-f) = 36KeV which is too small compared with 
the experimental value. Eq. ( |Tl| ) for c = 1 yields a decay width T{p — > n'y) = SlKeV 
in much better agreement with the data (see below). The experimental value of the pnir 
coupling or the KSRF relation [|T7| for the pTrvr coupling are used to calculate these widths. 

One can improve these equations by making the vector meson p unstable using the self- 
energy correction for the p propagator |T^ and using the KSRF relation |T^. This same 



result can also be obtained using the inverse amplitude for the vector form factor without 
assuming the KSRF relation. The p width obeys the KSRF relation as a consequence of 
the implementation of the unitarity relation [|1^]. The factor m^/ (m^ — s) is then replaced 
by a function Q{s) which is normalized to unity at s=0 and is deflned as follows: 

^^'^ ^ 1 - s/sr - 96^{(s - 4ml)H^^{s) + 2s/3} ^^^^ 

where f^, = 0.0936*61^, and sr is related to the p mass squared = 0.593GeV^ by 
requiring that the real part of the denominator of Eq. ( [TBI) vanishes at the p mass; H^^j^ls) 
is a well-known integral over the phase space factor: 




Ami , + V ^ - j s - Ami . . 2 



H^^is) = 2-2\\ -^In ^- + Z7r^/ s > 4m: 
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Ami — s I i n 

2-2W 1 arctan./ < s < Ami 

V s \ s — Ami 



s-Aml , JAml-s + ^/^ 
= 2-2\ ^In-^^ ;s<0 17 

Let us call the phase of Q{s) 6. Then Q{s) has the following phase representation: 

^/ N s 6(z)dz , , 

1] s = exp - / 18 

The phase 8 is exactly the elastic P-wave vrvr phase shifts as can be seen from Fig. 1. 
Alternatively one can use the experimental phase shift to calculate the function but 
the expression given above is most convenient. 

Other functions f2(s) normalised to unity at s = Sq can be expressed in terms of the 
function f2(s) by the simple relation Sq) = n(s)/n(so). 

as given by Eq. ( [T6| ) has a ghost pole at s = — 2.5.10^6*6^^ which is far away 
from the physical region relevant to our calculation and hence is irrelevant for our low 
energy calculation. 

The function defined here is the same as the inverse of the D-function given by 
reference |T0(] except for the definition of the p mass which is approximate there. 

In both approaches, the chiral symmetry limit should be defined as the limit of s, t, u 
tend to zero first and then — 0. This order should be respected because the branch 
point at s,t,u = Ami. ^Iso goes to zero in the chiral limit. Using this definition we could 
have calculated A in terms of A without using the large A''^ limit, but the difference is 
negligible as discussed previously. 

Replacing m^/ (m^ — s) by Q{s) in Eq. ([T0|) yields r(p 717) = A2KeV and with c=l 
(the hidden symmetry model with an additional assumption of a complete vector meson 
model for 7r° — > 77) in Eq. (|ll|) gives r(p — 717) = 95KeV. The difference between 
these values and those obtained previouly are just due to the p finite width correction. 
These results show the importance of the presence of the contact term. While the present 
experimental data on the r(p ^ iT'-f) are not settled, it is likeky that the result for the 
hidden symmetry model is favored (see below). 

With chiral symmetry broken, the pions acquired a finite but small mass, Eqs. (ffOl 



TT],IT2D become, respectively: 



and 



G^""^{s, t, u) = -Ms) + n{t) + n{u)} (19) 

3 

G^"'^"(s, t, u)2 = -^{n{s) + n{t) + n{u) - c} (20) 

o C 

G-'-{s,t,u) = ^{[fi(.)(l + ^;^)] + [s^t] + [s u]} (21) 
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Eqs. (p!9| , |20| , piD , do not satisfy, however, the elastic unitarity relation, i.e. the projected 
P-wave do not have the phase of the P-wave tttt interactions below 1 GeV JlQl as can 
be shown in Fig. 3. This result is not surprising because the multiple tttt scattering 
correction, which should be relevant for this problem, is not taken into account in these 
equations. The contact term model with c=l satisfies the phase theorem better than the 
pure VMD model due to the presence of the contact term which increases significantly 
the magnitude the resonance term compared with the background terms from the t and 
u channels. This leads us to a smaller correction using the following integral equation 
approach where the unitarity relation is explicitly built in. 



Eqs. (|T0|,|rT],|T2D can be considereds as the large Nc limit of the QCD. This is true 
because ~ Nc and the function Q{s) defined by Eq. ([T6|) becomes in this limit a simple 
pole. We shall elaborate this fact later in this article. 

In the following, we shall define the function G{s,t,u) at the symmetry point s = 
t = u = while the chiral symmetry limit of this function is the chiral anomaly A 
given by Eq. (^. How are they related to each other? There is no clean answer for this 
problem. Chiral Perturbation Theory could be used. The answer depends, however, on 



one parameter which is the scale parameter |rT|. We prefer to look at the large Nc 
limit to get their relation. 

Setting s = t = u = Om the chiral limit in Eqs.([TO|, |ll|), we have: 

_ 3 

G{s = ml,t = ml, u = ml)^X = X[l + -^ 1] (22) 

3 — c rrip 

This expression will be used in the following analysis. For c=l, we have A = 1.049A 
whereas the corresponding value for the one- loop ChPT ||^, assuming that the scale 
parameter /x^ = m^, is A = 1.053A which is insignificantly larger. In terms of A, with 
chiral symmetry broken but in the large Nc limit, one has: 

G.-=(..Va) = -^{!^V^ + !if^ + !2V^-c} (23) 

3 — c — s — t — u 



3.2 Vector Meson Model for Pion Form Factor 

The solution of the integral for the pion form factor with the assumption of the elastic 
form factor is 

Vis) = PnisMs) (24) 

where Q{s) is given by Eq.(|l^) and Pn{s) is a polynomial of degree n in s with real 
coefficients and Pn(0) = 1. 

For a given a set of the strong P-wave tttt phase shifts, the solution of the MO equation 
is not unique. One can multiply the solution Q{s) by a real polynomial to get a different 
set of solutions with different asymptotic conditions. The low energy constraint enables 
us to fix at least some coefficients of the polynomial. 

If one assumed P„(s) = 1, the square of the modulus of the pion form factor at the p 
mass is too small by about 30%, see Fig. 2, and the rms radius of the pion is too small 
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by 10%. Constraining the rms radius to be equal to its experimental value, we have to 
set [ig Pni-s) = 1 + 0.15(s/m2) or: 

V{s) = {l + 0.15^)Q{s) (25) 

The connection between this equation and the contact term was recently discussed 
The following integral equation for the process —>■ vrvr is more complicated because 
the integral equation involves both the right and left cuts on the real axis and hence the 
ambiguity of the solution is not simply the polynomial ambiguity and but is only related 
to it. It can only be obtained by solving numerically the integral equation as will be 
shown below. 



4 Integral Equation Approach using Elastic Unitar- 
ity Relation 

In this article, the process 77r — >■ tttt is studied using dispersion relation and elastic uni- 
tarity for the lowest partial wave. A similar integral equation of the type Muskhelishvilli- 
Omnes integral equation is obtained. The difference is that the integral equation to 



be treated here is much more complicated due to crossing symmetry; no exact solution 
has been found. We shall get the solution of this integral equation by an iterative proce- 
dure, but with the crucial property that the iterative solution for the lowest partial wave 



satisfies the phase theorem at every steps as required by unitarity |]I9|. As the solution 
of the M-0 equation is ambiguous by a polynomial, we find here a similar problem. But 
the ambiguity is not the same i.e. a new solution cannot be obtained by multipying the 
old solution by a polynomial. 

We start first by deriving the single variables s,t,u dispersion relation for 77r tttt; we 
then project out the P-wave amplitude where the rescattering effect is important due to 
presence of the low energy p resonance at 0.77 GeV. The rescattering effect is supposed 
to be negligible for higher partial waves because there are no resonances below 1.5 GeV 
for the two pions in F, H... waves. After solving the integral equation numerically, one 
should put the results obtained in to a form of the single variables dispersion relation. 
The crucial point is the single variable dispersion relation for the scattering amplitude do 
not satisfy the phase theorem, but its P-wave projection does. 

The integral equation can be derived using the technique of the Roy's equation for 
TTTT — s> TTTT Scattering |2^. We begin by writing a twice subtracted dispersion relation for 



G{s, t, u) at a fixed t. This dispersion relation can be shown to be valid in general. Using 
the same technique as that used in obtaining the Roy's equation, namely using the fixed 
t dispersion relation, crossing symmetry and keeping only the P-wave for the partial wave 
expansion of the absorptive part, one arrives at ( for an explicit demonstration of this 
equation see reference |1TT|| ): 



, ^ A . I £ „_„;^^,_,J + |. » «1 . |. « (26) 
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where the symmetry point in the problem is ai s = t = u = m^, and A is related to A as 
will be discussed later. The subtracted linear terms do not contribute because they are 
proportional to (s — m^) + (t — ml) + {u — m^) = 0. One can make a partial fraction of 
the dispersion integral to show that one can equally work well with the once subtracted 
dispersion relation which we shall use in the following: 

^/ , N r(s-"^^) ImGi(z)dz , r ,1 r T 
G{s,t,u) = A+ ^ ^ + + 

= A{s) + Ait) + Aiu) 

(27) 

where ImGi{z) is the imaginary part of the P-wave amplitude. 

The assumption of the dominance of the P-wave amplitude to get the integral equation 
can be experimentally checked by measuring the absence of the deviation from the sin^ 9 
angular distribution given by Eq. (|^). Should this assumption be incorrect, one could try 
to treat phenomenologically the contribution of the higher partial waves in the p region 
by some real amplitudes [|n[]. 

Let us call Gi{s) the lowest P partial wave projection of G{s, t, u) as given by Eq. (H). 
The elastic unitarity relation gives: 

ImGi (s) = Gi (s)e-^^(') sin 5{s) (28) 

where 5 is the P-wave vrvr phase shifts obtained from the available experimental data 
which show that they pass through 90° at the p mass as can be seen from Fig. 1. There 
is no measurable inelastic effect below 1.2 GeV. Projecting out the P-wave from Eq. (|27| ) 
and interchange the order of integration, we have: 

^ , ^ X s-ml f'^ Gi{z)e-'^^'''^ sm.5{z) , 3 /"^ ^ , , . , 

Gi s = A + ^ . ' ... ^-^dz + — Gi{z)e '^^'>smS{z) X 

TT J Ami {z - ml){z ~ s - le) 2-n JAml 

r 1 I z-a{s) + h{s) z - a{s) 4 

^b{s)^ b{s)^ ' ^ z - a{s) - b{s) ^ 6(s)2 3{z-ml)^ 

(29) 

The first derivative of Gi{s) at s = vanishes and its second derivatives with respect 
to s evaluated at s = is: 

d'Gijs) I _ 12 ImG,{z)dz ^^^^ 

The standard solution of the M-0 equation is ambiguous by a polynomial, but this 
problem is much more complicated here due to the symmetry of the s,t,u channels which 
leads to a much more complicated integral equation and hence it does not have the same 
type of ambiguity. Eq. (|29|) is a complicated integral equation. It is similar to, but more 
complicated than the Muskelishvili-Omnes (MO) type []12|, because the t and u channel 
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contributions are also expressed in terms of the unknown function Gi{s). It should be 
noticed that the first term has a cut from 4m^ to cxo and the second one has a cut from 
to — oo. For s > 4m^ the argument of the logarithm function in Eq. ( pQ]) never vanishes 
and hence this enables one to solve the integral equation by the following iteration scheme 
which converges very fast. 



4.1 Iterative Solutions 

As it was remarked above, the integral equation (^) has both right and left cuts. In 
setting up the iterative scheme, it is important to keep in mind that the final solution 
should be symmetric in the s,t,u variables as given by Eq. (pTf ). Because of this analytic 
structure, we can define an iteration procedure which consists in splitting Eq. (|29| ) into 
two separate equations: 

and 



3 TT J4:ml [z — m^){z — s — le) 



Tt'\s) = ^ + f r G«(.)e-^(^)sin5(.)x 

O ZTT J Ami 



r 1 I z-a{s)+b{s) z - a{s) 4 

^b{sy b{sy ' ^ z - a{s) - b{s) ^ b{sy 3{z-ml)^ 

(32) 

where i > I and Gi^ is the value of the function Gi(s) calculated at the z*'* step in the 
iteration procedure; the Born term T^^(s) is calculated at the i^'^ — 1 step. An iteration 
cycle is defined as a numerical calculation of both these two equations. 

The Born term is real for s > and has a left cut in s for s < 0. In writing Eqs. 
(^,^), care was taken to preserve the symmetry in the s,t,u variables for the function 
F{s,t,u) which requires us to spht symmetrically the subtraction constant A in Eq. (l29| ) 
into three equal pieces, one contributes to Eq. (^) the other two to Eq. 

The solution of the integral equation Eq. (|3TD is of the MO type [O: 



where: 



mz 



TT 



?>m.6{z)T^ '{z)dz 



Ami 



{z — ml){z 



s — te] 



(33) 

(34) 

i.e. this new function Q is normalized to unity at s = and Pn{s) are polynomial of 
nth degree with real coefficients. In the following, similar to Eq. (pSj), we only take the 
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first to terms in the polynomial and hence set: 



P„(s) = l + «^ (35) 

where a is a parameter which is related to the contact term c defined previously. The 
second derivative of the P-wave amplitude, defined by the sum rule Eq. (pOD, depends 
sensitively on the parameter a. 

Eq. ( P^D is not really typical solution written down for this type of integral equation. 
It is usually written in terms of the "driving" term Tg~^''(z). This procedure is not at all 
valid for the present situation, but we must modify it in order to get the final solution 



for the full amplitude with is completely symmetric in the s, t, u variables. Eq. (33) 
is written with this fact in mind. The first term on its RHS represents the VMD in the 
s-channel with or without contact term, the second term is the corresponding contribution 
from the t and u channels and the third term is the rescattering due to the final state 
interaction in the s-channel. 

At first sight one would think that the RHS of Eq. (|33|) does not have the P-wave 
phase 6. This is not so because we first note that the last integral can be separated into a 
principal part integral and a delta-function contribution which is purely imaginary. Then 
combining this delta-function contribution with the T^* ^■*(s) in Eq . (|33D , we have: 

Gf{s) = U{s,ml){^ + Tt'\s)Re\n"\s,ml)] + 

s-ml n^\z,ml)e'^^'hm6{z)T^^~^\z)dz 



TT Jiml {z-ml){z-s) 

(36) 

where P stands for the principal part integration. We have made the usual decomposition 
N/D for the partial wave amplitude, which can be shown to be quite general and inde- 
pendent of any dynamical scheme: e*'^sin5(s) = N{s)p{s)n{s,mD where N{s) contains 

only the left hand cut,D~^{s) = Q{s) only the right hand cut and p(s) = ^JJY--4m^Js) . 

Eq.(|36D shows indeed Gi\s) has the phase 6. 

One arbitrarily defines the convergence of the iteration scheme at the ith iteration 
step when | \ / \ Gi \ differs from 1 by less than 1% or so in the energy range 
from the two pion threshold to IGeV. (Alternatively one can also require that the ratio 
I I / I ^B^^"* I to be unity within an accuracy of 1%). 

Once the solution for the partial wave is obtained, one should return to the calculation 
of the full amplitude. This can be done by combining the T^~^^ Born term in Eq. ([33| ) 



with higher uncorrected partial waves (for rescattering) from the t and u channels to get 
the final solution: 

G«(.,t,«) = ^[{n(.,m^)(l + 3/(-^)(.))} + {{s ^ t)} + {{s ^ u)}] (37) 
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where the function I^^ denotes the multiple rescattering correction: 

TT J 4ml - ml) {z - s - ie) 

It is obvious that G^'^\s,t,u) do not have the phase 5 of the P-wave tttt scattering, but 
its P-wave projection does. This is so because, projecting out the I = 1 partial wave from 
Eq. (0), we arrive at Eq. (^) with T^~^(s) replaced by T^{s). Because of the assumed 
criteria for the convergence of the iteration scheme, T^~^(s) ^ T^{s), it is easily seen 
that Gi\s) has the phase 6, using the result of Eq. ( ^31) and Eq. (|36|) . The remaining 
higher partial waves / > 1 are all real because we have assumed that the strong final 
state interaction of the higher partial waves are negligible. The final solution Eq. (^) is 
completely symmetric in the s, t, u variables. 



5 Numerical Solutions 

We shall solve numerically our integral equation for various values of a defined by Eq. 
(p5|), corresponding to different values of the contact terms as discussed in Eq. (|ll|). We 
examine the following case: a = 0, 0.30, 0.50, 0.70. The iteration scheme can be done by 
first guessing a solution for Tb{s) corresponding to a chosen value of a. We can take a 
rather arbitrary first solution for this function. For example, for a given a we can take 
the t and u channel contribution to the Born term as: 

Tl{s) = rfcos^sin2^-{[^^tl^] + C''"^' ]} (39) 

which is independent of a. With this expression for the Born term, the iteration scheme 
can be started by calculating the solution of the integral equation as given by Eq. (PBf). 
The next step is to calculate the Born term by Eq. (|32D where there is no arbitrariness 
at this step. The iteration scheme is continued until the convergence of the solution is 
obtained. 

The number of the iterations depend on the original choice of the Born term i.e. 
how closed it is to the final solution. Even for a not so good approximation for T^^ as 
given by Eq. (|5UD, it is found that after one iteration one can already reach a reasonable 
approximation for the solution of the integral equation. 

Convergence to the final solution for a = 0.3, 0.5, 0.7, with a precision of the order of 
1% or better, could be achieved without using the iteration scheme if one chose a good 
expression for T^{s). For this purpose, one used instead the Born terms calculated from 
the following zero width contribution of the t and u channels of Eq. ( [TT| ) which depends 
on the strength of the contact term: 

^0. N 3 /-i ^ . 2^ A(l + a) ^^ m^-m^ ^ rn^ - ml 2a 

Tr s = - / rfcost^sm fc'^ \\ — '-^ + — \ 40 

^ 4 7-1 3 m^-t ^ ^rn^-u^ I + ^ ' 



11 



where we have used the large Nc relation: 

"=3^ (41) 

We only make use of this relation here to calculate the Born terms and also to get the 
relation between A and A. 

This relation is obtained in the limit of of a narrow p width. As it was discussed above, 
this limit is obtained when we let the number of color oo. In this limit, as it will be 

shown below, Eq. (^) becomes Eq. (0) or Eq. ([T^). The pure VMD model corresponds 
to a = and for models with the contact term c=l, e.g. the hidden symmetry model 
0, a = 1/2. (More precisely, the hidden symmetry model with no contact term in the 
pseudoscalar mesons decaying into 2 gammas requires c=l). For the real situation where 
Nc = 3 there should be substantial correction to this relation. 

For a = 0, without an iteration of the integral equation, one can only get a precision 
of the order of 5%. Similarly for other values of a, using Eq. (|D) for the Born terms, 
and without going through the iteration scheme, one can already achieve a precision of 
better than 1% using only Eq. (PT|). 

For a = 0, after 5 iterations a precision of better than 1% is reached. For other values 
of a= 0.3, 0.5, 0.7 a precision of better than 10~^% is reached after only 4 iterations. 
These numbers indicate the rescattering effect is much more important for a = and is 
much less important for other values of a. 

The slow convergence of the iteration scheme for a = is due to a large violation of 
the phase theorem at the zeroth order as shown in Fig. 3. For other cases the violation 
of the phase theorem is not so serious and even without the iteration scheme one can 
already get the approximate solution accurate to better than 1% by solving directly the 
integral equation as discussed above. 

Instead of parametrizing our solution by the value of a, it is more physical to describe 
the solution as a function of the width r(p TC'-f). This quantity is not unambiguous and 
will be defined in the following section. It is denoted by r(p ^ n^y) using our definition 
while the corresponding partial width using the usual Breit-Wigner parametrization is 
denoted by r(p — > 'n''j)bw with the value of the p mass the same as our defintion, i.e 

= 0.593GeV^ and r(p — > tctt) = 0.156GeV. There is a substantial difference for the 
values of r(p nj) using these two definitions. In Table 1. r(p — * nj), T{p ^ T^l)hw 
and the second derivative of the P-wave amplitude at s = m\ are given as a function of 
the numerical value of a. 

Corresponding to various values of a the square of the absolute values of the P-wave 
amplitudes Gi(s) are plotted against the energy square s, in the unit of A, are given in 
Fig. 4 and Fig. 5. It is seen that the larger the values of a the higher is the maximum 
values of the P-wave amplitudes. At lower energy, the absolute value of the magnitude of 
the P-wave amplitude also increases with a until a = 0.5 and then decreases with larger 
values as shown in Fig. 5. 

In Fig. 6 the modulus of the ratio 3Gi{s)/Q{s) is plotted against the energy squared 
s(GeV^); this ratio indicates the deviation from the Breit-Wigner form as given by the 
function fl{s). 
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For various values of a, accurate values (to better than 1%) of the modulus of the 
P-wave amplitude from the two pion threshold to 1 GeV can be obtained by using the 
modulus of f2(s) divide the function C(s) given in Table 2: 

G,{s)=U{s,ml)C{s) (42) 

We are also interested in finding out what are the corrections due to the multiscattering 
effect in the VMD approximation to the function A(s) defined by Eq. (|2^) and given by 
Eq. O: 

Ms) = '-^Jis) (43) 

or 

As) = ^il + 3r-\s)) (44) 

for a value of i attained at the end of the iteration of the integral equation. In Fig. 7 the 
modulus of J(s) is plotted against the energy squared s in the unit of GeV^. If there were 
no corrections to the VMD model, J(s) would be unity. It is seen that the corrections are 
most important for the case of a = 0. 

In Fig. 8 the phase of A(s) is also plotted against the energy squared s for various 
values of a and compared with the P-wave tttt phase shifts. 

In Fig. 9, I G{s, t, u) p with cos 6* = and cos 6* = 0.75 are plotted against s (in GeV^) 
are shown for the special case a = 0.5. Figs, with other values of a and cos 6' are not 
shown because they are quite similar to the Fig. 9. Therefore the higher partial waves 
are completely negligible for an energy below 1 GeV. 



6 Comparison with experimental data and other the- 
oretical works 

Our calculation can be compared with the experimental data at low and high energy. At 
low energy, the only experimental data available is given by reference [Q. From Fig. 5, at 
an energy s ~ O.lQGeV^, corresponding to the average energy measured in the reference 
1^ we have, for a=0.5, | Gi(s) |~| G{s,t = u) |= 1.15A which is about one standard 
deviation smaller than the measured value (1.29 ± .09 ± .05)A. It is clearly important to 
improve the precision of this experiment. 

At higher energy, the experimental cross section for 77r — >■ tttt are usually analyzed in 
terms of the Breit Wigner formula: 

da_ ^ 2ATIS mlT{p ^ 27r)r(p ^ 717) ^^^^ 
ds {s — (m^ — s)2 + m^r^(s) 

This formula is usually not accurate because it either neglects the contribution of the 
part of the amplitude from the t and u channels or assumes that the cross section can be 
fitted with a Breit Wigner form which may not be true. Furthermore the maximum of 
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the modulus of the P-wave amphtude is shifted significantly toward lower energy which 
complicates the analysis of the experimental data using Eq. (|45|) . 

The result of our calculation shows that, at the maximum of the absolute value of the 
P-wave amplitude, the phase of the amplitude is not 90°. The only reasonable method 
which we can find acceptable is to define the p mass as the value of s when the phase 
of the function which is the same as the experimental P-wave vrvr phase shift passing 
through 90°. Its width is proportional to the inverse of of the derivative with respect to 
s of the cot 5 aX s = rn?^: 

= cot 5{s) |^=„2 (46) 
mpL p as 

With this definition the p width as given by Eg . ([I6|) is 0.156 GeV. One then could use the 
P-wave cross section at s = to calculate the r(p —>■ 717) width using Eq. (|45|). The 
values obtained is denoted by T{p 717) and would approximately be 10% lower than the 
value obtained from using the maximum observed 77? — vrvr cross section in combination 
with the Breit-Wigner formula, Eq. (|45|) which is now denoted by T{p 7r7)b^ (see Table 

!)• 

Using our method, we could even integrating the measured cross section on either side 
of the p mass by 0.1 GeV in order to improve the experimental accuracy without changing 
its value by more than 1%. This precision would not be possible if one did the calculation 
with the p mass as the value of the maximum cross section. 

The present experimental results are not consistent with each other. The more recent 
published experimental results by Caparo et. al. gave the value for T{p —>■ 717) = 81±4±4 
KeV 1^^, wheras earlier results by Huston et. al. gave a lower value These two 



experiments are the Primakoff-like experiments using a high energy charge pion beam on 
a heavy target. The experimental result from e~^e~ reaction gives a higher value for the 
p Tf-f width but experimental results contain large error. 

A more recent unpublished result using photo-production of a pair of pions off the 
nucleon target yields T{p — > itj)},^} = 96 ± 12KeV [^]. Unlike two previous Primakoff 
experiments, this experiment might have some difficulties with isolating the data corre- 
sponding to the one pion exchanged diagram from the background effect; one must also 
take into account of the fact that the exchanged pion is off its mass shell. 

Due to the lack of experimental informations on the second derivative of the P-wave 
amplitude at s=0 or the parameter a (see the Table 1), we cannot predict the solution of 
the integral equation to get the r(p — > ti'j) width. 

Corresponding to a naive pure VMD model without a contact term, from Table 1, 
our calculation with a = yields a width of T{p — >• nj) = 50.7KeV or T{p 7i'j)bw = 
57.8KeV, whereas corresponding to the hidden symmetry model with c=l, our calculation 
with a = 0.5 yields T{p 717) = 8A.3KeV or T{p 717)^^ = 91.8KeV. 

With a = 0.5 the value for r(p TT'~f)bw is somewhat smaller than the value of 96 KeV 
obtained by Hannah using the Pade or the inverse amplitude methods for ChPT two loop 
amplitude which was calculated numerically For this special value of a, one would 
also recover the main result of the hidden symmetry model with c=l. The low energy 
parameters C and D, corresponding to the first and second derivatives of the function 
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A(s), Eq. (^3]), defined and evaluated by Hannah |TT[ are in agreement to an accuracy of 



2% witli tliose from tlie our integral equation approach. 

The difference between this work with that of Hannah is presumably due to the in- 



terpretation of the Eq. (45), the treatment of the multiple scattering effects and also 
to the interpretation of the contact term. Hannah's work shows the importance of the 
resummation of the perturbation series by the inverse amplitude or the Pade methods. 

There is a similar treatment of this problem by Holstein. Holstein solution was ob- 
tained by taking the product of the three functions: 

Gf (s, t, u) = AP„(s, t, u)n{s)n{t)n{u) (47) 

where Pn{s,t,u) is a polynomial in s, t, u constructed in such a way that this equation 
has the same low energy limit as that given by ChPT. The merit of this equation is the 
phase theorem is explicitly obeyed as can be seen by projecting out the P wave from this 
equation. But this equation is not right because all higher partial waves such as F, H 
...all have the p resonance or they all have the phase of the P-wave phase shifts which 
are not correct. The singularity associated with the multiple scattering effects which are 



present in our integral equation approach do not contain in Eq. (|47|) . All possible solutions 
which can be written in terms of the product form of three functions in s,t,u variables will 
have this problem. One exception is the problem involving three hadrons with two light 
particles having no interactions between them but they interact with an infinitely heavy 
target. 

Holstein's solution yields a comparable value for T{p 'n'~^)bw as our solution with 
a = 0.5. This result is not surprising because the value of the second derivative of his 
solution at s = is also comparable with ours. His solution can be fixed by projecting 
out the P-wave imaginary part and put it in Eq. (^61) to provide the necessary corrections. 

It should be reemphasized that our result is not in the product form as in Eq. p7|) but 
as a sum of the three identical functions with the interchange of s,t,u variables Eq. (|26|) 
and Eq. (|37|) . It is a direct consequence of the fixed t dispersion relation, using crossing 
symmetry and neglecting the contribution form the higher partial waves at low energy in 
the absorptive parts. 



7 Importance of the Multiple Scattering Correction 
as a Function of the p Width 

Our formulation of the problem 77r tttt is quite useful in understanding the importance 
of the multiple scattering effect as a function of the p width. We have previouly stated that 
in the large limit the multiple scattering effect should vanish and we should recover the 
VMD models as given by Eqs. (|TD|JTT|,|T2|) . In order to see that this statement is correct, 
it is sufficient to study the correction factor J(s) defined by Eq. (|43|) as a function of the 
p width. It is sufficient to study this question for the case when a = 0. In Fig. 10, the 
modulus of the function J(s) is plotted against the energy squared when the p widths are 
increased or decreased by a factor of 4 when f^^ is changed by a factor of 2. This can be 
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seen from examining the definition of the function of Eq. ([T6|). The result is the 

multiple scattering effect increases as the p width increases, and decreases as the the p 
width decreases. It is easy to see that for a zero width resonance, the correction factor 
J(s) is unity. 

The result of this section provides also some arguments for the neglect of the multiple 
scattering effect in the study of < 7 | Svr > and < Svr | 7r7 > with the Svr resonating as 
the uj state because of the extreme small width of this resonance [ pT[] . 

8 Conclusion and Acknowledgment 

We have studied in this article the scattering of 77? tttt using the integral equation 
approach. Because the second derivative of the P-wave amplitude at s = is not 
known, the maximum cross section for this process cannot be predicted with reliability. 
The solution of the integral equation is ambiguous and depends on the second derivative 
of the P-wave amplitude at s = m^. This problem is similar to problem of the contact 
term in the usual VMD model. 

If the ambiguity of the solution of the integral equation could be interpreted as the 
imperfectness of the elastic unitarity relation in describing the low energy phenomena, 
then one would have to be satisfied with a precision of the order of 15% in amplitude for 
the pion form factor calculation [Q. This inadequacy could then be removed using the 
knowledge of the pion r.m.s. radius [^. 

For the 77? tttt calculation, this inadequacy could become more serious due to 
the existence of singularities associated with the t and u channels. Furthermore, the 
corresponding first derivative of the P-wave amplitude vanishes due to the symmetry of 
the problem, and hence we could only use the knowledge of the second derivative at 
s = to fix up the inadequacy of the elastic unitarity relation. This last parameter 
could not be precisely measured and hence we could not predict with certainty T{p ^ tij). 

We show in this article that there is a one to one correspondance between the contact 
term in the VMD model and the ambiguity associated with the parameter a in our integral 
equation approach. This parameter plays the role of the second derivative of the P-wave 
amplitude at s = m^. 

It is a pleasure to thank Torben Hannah for sending me a copy of his paper before 
publication and for useful discussions. 
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a 


r(p — >• 7r7) KeV 


r(p ^ 7r7)fe^ KeV 


\~'(PGi{s)/ds''{s = ml)GeV-^ 


0.00 


50.7 


57.8 


4.88 


0.30 


68.0 


76.4 


5.30 


0.50 


84.3 


91.8 


5.64 


0.70 


103 


111 


6.00 



Table 1: Solution of the P-wave amplitude for the 'yir —>■ vrvr integral equation as a 
function of the parameter a. The second column is the r(p ^ n'-f) in KeV according to 
the definition given in the text. The third column is the r(p 'n''y)bw in KeV using the 

Breit Wigner formula Eq. p5|) evaluated at the maximum of the cross section; The 
fourth column is the second derivatives of the P-wave amplitude at the point s = in 

the GeV'"^ unit. 



a 


C{s) 


0.00 


2.457s3 - 1.662^2 + 1.162s + 0.2813 


0.30 


0.808^2 + 0.667s + 0.3250 


0.50 


-0.759s3 + 1.278s2 + 0.705s + 0.3267 


0.70 


-1.015s3 + 1.183s2 + 0.898s + 0.318 



Table 2: Relation between the P-wave amplitude Gi(s) and f2(s,m^) as given by the 
function C(s) defined by Eq. (^) given in text. 
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Figure Captions 

Fig.l : The phases of the function Q[s) in degrees (vertical axis) are given as a function 
of s{GeV'^). The experimental data are taken from references p8| , p9| , pO|] . 

Fig. 2 : The square of the modulus of the function Q{s) (vertical axis) is given as a 
function of s in GeV^ (dashed line). The square of the modulus of the pion form factor 
V{s) as given by Eq. (pSD (solid line). Experimental data are taken from references 

mm- 

Fig. 3 : The P-wave strong vrvr phase shifts (vertical axis) as a function of the energy 
are shown by the solid line. The projected P-wave amplitude phases of the 77r° — >■ tt+tt" 
as given by the VMD model without contact term, Eq. (0), are given by the long dashed 
line; with the contact term c=l, Eq. (|20|) , are given by the short dashed line; with the 
contact term c=l, Eq. (pi]), are given by the medium dashed line. 

Fig. 4 : Plot of the square of the absolute value of the P-wave amplitude | Gi{s) p in 
the unit of A as a function of s in GeV^ for a = 0.5 (solid line); a = 0.0 (short dashed 
line); a = 0.3 (medium dashed line); a = 0.7 (long dashed line). 

Fig.5 : Same as Fig. 5 but with O.OSGeV^ <s< 0.25GeV^. 

Fig. 6 : Plot of the ratio 3Gi{s)/Q{s,m'^) in the unit of A for various values of a. The 
curves are same as in Fig.4. 

Fig. 7 : Plot of the absolute value of the function J(s) defined by Eq. ( ^3] ) in the unit 
of A as a function of s in GeV^. The curves are same as in Fig.4. 

Fig. 8 : Plot of the phases of the function A(s) defined by Eq. (^) as a function of 
s in GeV^. The solid line represents the P-wave strong vrvr phase shifts; the dotted line, 
a = 0.5; the short dashed line, a = 0.0; the medium dashed, a = 0.3, the long dashed 
line, a = 0.7. 

2 

Fig. 9 : Square of the modulus of G(s,t,u) in the unit of A with cos6' = (long dashed 
line), with cos 6' = 0.7 (short dashed line) and the square of the modulus of the P-wave 
amplitude Gi{s) (solid line) for a = 0.5 as a function of s in GeV^. 

Fig. 10. : Plot of the absolute value of J(s) (vertical axis) defined by Eq. (|i3| ) in the 
unit of A vs s (GeV^) for a = for various values of the p widths. The solid curve 
represents Pp = 0.156GeV, the short dashed curve, Pp = 0.0396*61^, the long dashed 
curve, Tp = 0.624Gel^. 
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